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Structure and lattice dynamics of the (001) surface of S^CuC^Ck have been studied by helium 
atom scattering (HAS). Analysis of diffraction patterns obtained by elastic HAS revealed a surface 
periodicity consistent with bulk termination, and confirms that the surface is non-polar and stable 
which favors a SrCl surface termination. Bulk and surface lattice dynamical calculations based on 
the shell-model were carried out to characterize the experimental phonon dispersions obtained by 
inelastic HAS. No experimental surface mode was observed above 200 cm -1 . Comparison between 
the experimental data and theoretical results for two different slabs with SrCl and Cu02 terminations 
showed that the experimental data conforms exclusively with the SrCl surface modes. 

PACS numbers: 



INTRODUCTION 



Surface lattice dynamics of metals, semiconductors, 
simple insulators and layered structures^ have been stud- 
ied extensively . However, high temperature super- 
conductors (HTSC), their parent compounds and re- 
lated lamellar cuprates are more complex in structure 
and dynamics and the study of their surface dynami- 
cal properties has been limited to only a few cases 2 ^. 
Because of the established relation between phonons 
and superconductivity in conventional superconductors, 
bulk phonons of the HTSC and their parent com- 
pounds received considerable attention in the early years 
following their discovery. Optical spectroscopies (in- 
frared (IR) and Raman anc j neutron scattering 
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were used to study the lattice dynam- 



ics of these compounds. Although it is believed that 
spin interactions play an important role in the super- 
conducting properties of HTSC cuprates, the role of the 
electron-phonon coupling still remains a puzzle. For ex- 
ample, the isotope effect in the optimally doped cuprates 
is very smalUi, but away from optimal doping a large 
isotope effect has been observed^. In addition, the 
highest longitudinal optical (LO) phonon branch soft- 
ens and shows line-width broadening when the parent 
compounds of several HTSC, such as La2- a; Sr a ;Cu0 4 15 i 16 
and YBa2Cu306+j^£, are hole-doped to a superconduct- 
ing state. Furthermore, in recent angle resolved photoc- 
mission spectroscopy (ARPES) experiments on several 
HTSC materials, an abrupt change of electron velocity at 
50-80 meV has been observed^ which is interpreted to 
be caused by the coupling of electrons to the LO phonon 
modes in the same energy range. Whether these effects 
suggest that electron-phonon coupling in HTSC can drive 
electron-pairing or it is just a secondary mechanism that 
occurs naturally in an interacting system without affect- 
ing the superconducting pairing mechanism is still open 
to discussion and interpretation 19 . Therefore the role of 



phonons in HTSC is still not completely understood and 
electron-phonon coupling can still effect the supercon- 
ducting state. For this reason lattice dynamical studies 
of these compounds and related cuprates are still impor- 
tant. 

In this paper we report on the study of the structure 
and dynamics of the (001) surface of Sr2Cu02Cl2- In 
the bulk, this compound is a layered perovskite with 
the body centered tetragonal structure of K2MF4 type 
and space group I4/mmm(DlJ l )2^. It consists of Cu02 
sheets separated by double layers of SrCl. Neutron 
diffraction studies^ , show that the lattice parameters for 
Sr2Cu02Ci2 at room temperature are a = 3.9716 A and 
c = 15.6126 A. The crystal structure of Sr2Cu02Cl2 is 
shown in Figure Sr 2 Cu02Ci2 is also isostructural to 
the high-temperature phase of La 2 Cu04, but with the 
out of plane oxygen ions at the apices of the CuOe octa- 
hedra replaced by CI and the La by Sr. 

However two clear differences between Sr2Cu02Cl2 
and La2Cu04 stand out: unlike La2Cu04, which un- 
dergoes a structural phase transition from tetragonal 
to orthorhombic symmetry at about 530 K 21 involv- 
ing rotations of the CuOe octahedra, no evidence of a 
structural distortion has been observed in Sr2Cu02Cl2 
down to 10 K22i22. Also in contrast to La2Cu04, which 
can be doped by substitution of La with a divalent 
ion like Sr +2 and shows a superconducting behavior 
below T c ~ 40K, neither electron nor hole doping of 
Sr 2 Cu02Cl2 has been possible^, although the similarly 
structured Ca2Cu02Ci2 can be doped successfully by Na 
with T c ~ 28 and very recently superconductivity 
at T c ~ 30 K has been reported in apical oxygen doped 
Sr 2 Cu0 2+ ^Cl2-^. 

Neutron diffraction experimentsS2iSSiSL22iS£ have been 
used to explore the similarities and differences between 
the magnetic properties of Sr2Cu02Cl2 and La2Cu04 
and to understand the magnetic interactions in these 
two-dimensional (2D) spin-1/2 systems. In particular, 
it is shown that both Sr2Cu02Cl2 and La2Cu04 are 
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FIG. 1: Crystal structure and positions of the atoms in 
Sr2Cu02Cl2. Atomic positions are taken from Ref. 



antifcrromagnetic (AFM) insulators with Neel temper- 
atures (Tjy) of about 250 K and 300 K, respectively, 
and have similar intralayer exchange coupling J of about 
130 meV; in the paramagnetic phase both are good ex- 
amples of 2D square lattice quantum Heisenberg anti- 
ferromagnets (2DSLQHA)SiiS£. The main difference be- 
tween the magnetic properties of these two compounds is 
that the orthorhombic distortion in La2CuC>4 allows for 
a Dzyaloshinsky-Moriya (DM) antisymmetric exchange 
interaction in the spin Hamiltonian which is absent in 
Sr 2 Cu02Cl2 because of the stability of its tetragonal 
structure to the lowest temperatures measured- 1 . The 
DM interaction in La2Cu04 leads to a weak ferromag- 
netic moment perpendicular to CuC>2 planes^. How- 
ever, the weak inter-layer coupling causes the moments 
of successive layers to order antiferromagnetically, so the 
weak ferromagnetism is hidden in the Neel state at zero 
fielcM.. 

In light of the above discussions, and the fact that 
Sr2Cu02Ci2 has the advantage that it can be eas- 
ily cleaved, a study of the surface lattice dynamics of 
Sr2Cu02Cl2 can not only provide information about the 
bulk phonons of this compound, which are currently lim- 
ited to a few optical spectroscopy studies, but also pro- 
vide insight as to the modifications effected at the surface 
to the vibrational modes of the bulk for La2Cu04, where 
the bulk dynamics has already been studied but no sur- 
face phonon study exists. 



II. EXPERIMENTAL PROCEDURE 

The single crystalline samples of Sr 2 Cu0 2 Cl2 used in 
our experiments were grown at MIT crystal growth fa- 
cility. The samples were plate like with the c-axis per- 
pendicular to the sheets of the plate. They were cut 
to dimensions of about 7x7x0.5 mm and mounted on 
a copper sample holder by silver conducting epoxy. A 
cleaving post was then attached to the sample surface, 
also by the epoxy. The prepared sample was mounted on 
a manipulator which accommodated xyz motions, as well 
as polar and azimuthal rotations. The in situ ultra-high 
vacuum (UHV) cleaving of the sample was implemented 
by the top-post method by knocking off the post on the 
surface. Immediately after cleaving, the quality of the 
long range ordering on the surface was checked by low 
energy electron diffraction (LEED). The LEED images 
were sharp spots and ordered in a square pattern, with 
no discernable satellites. 

Elastic and inelastic He-atom scattering (HAS) tech- 
niques were used for measuring the structure and dy- 
namics of the surface. Since about 1975, novel surface 
spectroscopic methods based on the interaction of ther- 
mal energy neutral helium atoms with solid surfaces have 
become possible due to the progress in combining high- 
resolution He beam production systems with UHV tech- 
niques. Because of their low kinetic energies (10 - 100 
meV), inert He- atom beams probe the topmost layer of 
the surface in a non-destructive manner. In addition, 
momenta of thermal-energy He atoms are of the order 
of the surface phonon momenta, making HAS suitable 
for scanning the entire surface Brillouin zone. Therefore 
helium atom scattering is our method of choice for the 
study of the structure and vibrations of the surface. 

A mono-energetic supersonic beam of thermal-energy 
neutral He atoms is generated by a nozzle-skimmer as- 
sembly and directed to the scattering chamber through 
various collimating slits. The energy of the He beam can 
be varied from Ei = 25 meV to Ei = 65 meV by changing 
the nozzle temperature from 110 K to 300 K. The tem- 
perature is controlled and measured by a diode sensor 
together with a temperature controller (Scientific Instru- 
ments, Model 9600-1). Polar rotation of the sample al- 
lows the variation of the He beam incident angle 9i with 
respect to the surface normal. The scattered He beam 
is collected by an angle resolved detector mounted on a 
two-axis goniometer which allows the scattering angle 9f 
to be varied independently. The detector is comprised of 
an electron gun and a multichannel plate (MCP) electron 
multiplier. The electron gun generates a well-collimated, 
monoenergetic electron beam crossing the He beam at 
right angles. The energy of the electron beam is tuned 
to excite the He atoms to their first excited metastablc 
state (2 3 S He*) upon impact. Deexcitation of a He* atom 
at the surface of the MCP leads to the ejection of an elec- 
tron which generates an electron cascade that is then col- 
lected by the anode of the multiplier. By electronically 
pulsing the electron gun a gate function is created for 
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time-of-flight (TOF) measurements in the inelastic HAS 
mode. The details of the detection scheme are given in 
Ref|23- All the measurements are performed with the 
sample surface at room temperature. 

By writing the He-atom wave vector as k = (K,k z ), 
where K is the component parallel to the surface, con- 
servation of momentum and energy for in-plane He scat- 
tering can be expressed as 



AK = k f sintf/ - k t sin 6^ = G + Q, 



and 



AE = Hlo(Q) = E f - Ei 



2m 



(1) 



(2) 



where subscripts i and / denote incident and scattered 
beams, respectively and AK is the momentum transfer 
parallel to the surface. G is a surface reciprocal lattice 
vector, Q is the surface phonon wave vector and hto(Q) is 
the corrsponding surface phonon energy. By eliminating 
kf from the above equations one obtains the so-called 
scan curve relations which are the locus of all the allowed 
AE and AK as dictated by the conservation relations: 



AE = E, 



sin Oi + AK/k t 



- 1 



(3) 



The intersections of these scan curves with the phonon 
dispersion curves define the kinematically allowed inelas- 
tic events for a fixed geometric arrangement. Thus, by 
systematically changing Ei, Oi and Of the entire disper- 
sion curves can be constructed. 
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FIG. 2: Typical diffraction patterns afong (100) and (110) 
azimuth directions for cleaved Sr2CuO2Cl2(001) surface. The 
intensity scale is logarithmic. 



III. RESULTS 

A. Elastic He Scattering 

Figure[21shows two typical He diffraction patterns from 
the (001) surface of Sr 2 Cu0 2 Cl 2 along (100) and (110) 
symmetry directions for 0i — 40° and 35°, respectively. 
The curves are shown as a function of the parallel mo- 
mentum transfer AK. The surface periodicity is ob- 
tained from the positions of the diffraction peaks along 
each direction resulting from the Bragg relation 

AK = nG {W0) , or AK = nG {no) , (4) 

where integer n indicates the order of diffraction. From 
these and the position of Bragg peaks obtained from 
other diffraction patterns taken at various incident an- 
gles, the value of a = 3.97 A for the surface lattice pa- 
rameter was obtained. This value is very close to the 
lattice parameter reported for the bulk 20 . This, and the 
fact that no systematic satellite peaks in the diffraction 
patterns were observed, suggest a bulk like termination 
at the surface and indicates surface stability against re- 
construction. 



In general, cleaving the Sr 2 Cu0 2 Cl 2 can result in two 
different surfaces with different compositions as shown 
in Figure 01 SrCl termination of the crystal results in a 
non-polar surface which, according to Tasker's classifica- 
tion of ionic surfaces^, consists of a neutral repeat unit 
of SrCl - Cu0 2 - SrCl, and is stable. However, termina- 
tion at a Cu0 2 layer, results in a repeat unit of Cu0 2 - 
SrCl which has a non-vanishing dipole moment perpen- 
dicular to the surface. Such a polar surface cannot be 
stable because of the divergence of the electrostatic en- 
ergy contribution to the surface energy. Hence, such ter- 
mination would be allowed only if the surface undergoes 
a reconstruction resulting in cancellation of the dipole 
moment perpendicular to the surface. The apparent sur- 
face stability against reconstruction on the (001) surface 
of Sr 2 Cu0 2 Cl 2 suggests that the top-most layer is indeed 
SrCl and not Cu0 2 . This is also supported by previous 
LEED and x-ray photoemission spectroscopy studies of 
the Sr 2 CuO 2 Cl 2 (001) surfac e) 34 ! 35 ' 36 which suggest that 
SrCl is the top-most layer. 

From the intensity analysis of the diffraction peaks, 
one can obtain the surface corrugation function. A two 
dimensional surface corrugation consistent with the 4mm 
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symmetry of the surface was chosen as 
/ 2tt(x + y) 



{2irx\ /27nA 

cos + cos 

V a J \ a J 



+ C2 



cos 



COS 



2w(x - y) 



(5) 



where Ci and C2 denote the corrugation amplitudes and 
are taken as model parameters. Values of £1 = 0.01 A 
and C2 = 0.15 A were obtained by fitting Bragg peak in- 
tensities calculated within the framework of the eikonal 
approximation^! and using an attractive potential well 
depth of D — 10 meV to the experimental intensities. To 
back up this calculation we also obtained the corruga- 
tion amplitudes by calculating the surface charge density 
p{x,y,z), since the corrugation function is the locus of 
the classical turning points for He atoms at a fixed sur- 
face charge density contour. p(x,y,z) was constructed 
as a superposition of the individual atomic charge den- 
sities at a SrCl surface. The atomic charge densities for 
Sr and CI were derived using the atomic wave functions 
of Herman and Skillman— . The values of £1 and £2 cal- 
culated at a constant contour of p(x,y,z) at a turning 
point of 2.7 A above the surface were 0.01 A and 0.147 
A, which are very close to the values obtained by the 
eikonal approximation. 



B. Surface Phonons 

Surface phonon energies were measured along the A 
and £ directions of the surface Brillouin zone (SBZ), 
shown in Figure Two typical TOF spectra showing 
diffusive elastic, as well as inelastic peaks are shown in 
Figure The data were collected using beam energies 
in the range of 25 - 65 meV. The high energy He beam 
(Ei ~ 65 meV) was used to probe the high-frequency 
surface phonon modes. The energy resolution for the 
lowest energy beam was about 3 meV. Phonon energies 
and momenta were calculated from the TOF peak po- 
sitions, after a Gaussian fit to the peaks, using Eqs.QJ 



and J2J. To characterize the ensuing phonon dispersion 
points lattice dynamical analysis for the bulk and surface 
(slab calculation) were carried out as outlined in the next 
section. 



IV. LATTICE DYNAMICAL ANALYSIS 
A. Static Equilibrium and Bulk Lattice Dynamics 

In lattice dynamical studies of the cuprates and HTSC 
parent compounds, the construction of a dynamical ma- 
trix based on pair potentials rather than mere force con- 
stants is far more advantageous in the sense that the 
underlying physics is transparent and many of the po- 
tential parameters for similar pair interactions can be 
transferred from one compound to another^S^. This 
is specially useful for compounds like S^CuO^C^ for 
which experimental phonon frequencies are limited to 
four infrared (IR)-active phonon modes at the center of 
the Brillouin zone (BZ)2ii22, and no neutron scattering 
data for phonons are available. However, experimental 
and theoretical studies of bulk lattice dynamics of the 
isostructural La2Cu04 are extensive and contain detailed 
information about relevant model potentials^Siiiii 2 ^. 
Therefore, for some of the pair interactions which are 
similar in both compounds (such as Cu-0 and O-O), the 
existing potential parameters for La2Cu04 can be trans- 
ferred to Sr 2 Cu02Cl2. In Sr 2 Cu02Cl2, with seven atoms 
per primitive cell, more model parameters can be de- 
termined by fulfilling the static equilibrium conditions 
which state that the forces on the particles in their equi- 
librium positions should vanish. This treatment also en- 
sures the consistency of the static and dynamical prop- 
erties of the crystal. It should be noted, however, that 
satisfying static equilibrium conditions, does not guaran- 
tee any dynamical stability, namely ensuring the reality 
of the phonon frequencies throughout the BZ. 

The model used for both static equilibrium cal- 
culations and lattice dynamics (shell-model) of bulk 
Sr2Cu02Cl2, incorporates two-body central potentials, 
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FIG. 5: Two typical TOF spectra along (100) and (110) direc- 
tions (solid lines). Dashed and dotted lines are the Gaussian 
fits to the peaks. 



namely Coulomb potential V§(r) = ZiZje 2 /r for the 
long-range interactions and either Born-Mayer (BM) 
Vg M (r) = a,j e~ b v r or Buckingham V?(r) = ay e~ b ^ r - 
Cij/r 6 potentials for the short-range interactions. Short- 
range interactions were limited to nearest neighbor pairs 
for Cu-O, O-O, Sr-O, Cl-Cl, Cu-Cl, and O-Cl, and to 
nearest and next nearest neighbors for Sr-Cl. Potential 
parameters were obtained from existing literature for Cu- 
O, O-O, Sr-0 and Cl-Cl. Static equilibrium conditions 
provided four equations and were based on the methods 



TABLE I: Potential parameters for pair interactions in 
Sr2Cu02Cl2 and shell- model parameters. 







Potential Parameters 




Ions 


a (eV) 


U \r\. J C 


^ev a j 


Cu-O a 5814.375 


4. f OZ 


U 


o-o a 


1146.25 


O.Z ( y 


U 


Sr-0 


1950 


2.978 6 





Cl-CT 18498.3 


3.65 


99.87 


Sr-Cl 


62037.2 


4.064 





Cu-Cl 103.55 


3.0 





O-Cl 


1161.75 


3.0 









Shell-Model Parameters 




Ion 


Z(e) 


Y (e) K (eV/ A 2 ) 


Cu 


2 


2.6 


40 


Sr 


2 


1.9 


30 


O 


-2 


-3.1 


24 


CI 


-2 


-1.8 


5 


a Ref. 
6 Ref. 
c Ref. 


45 
46 



TABLE II: Comparison of the four measured infrared-active 
frequencies for Sr2Cu02Cl2 at the center of the Brillouin zone 
(q — 0) with the values calculated by the shell model in this 
work. All the frequencies are in units of cm -1 . 



Measured Frequencies 



Mode 


Ref. 41 


Ref. 42 


Calculated Frequencies 


E U {1) 


525 


512 


522.41 


E u {2) 


351 


339 


334.37 


E u (3) 


176 


173 


182.67 


E u (4) 


140 


138 


143.51 



developed in Refs. 03 andEfl. 

Table [I] lists all potential parameters obtained either 
from the literature or by satisfying equilibrium condi- 
tions. These parameters were then incorporated into a 
shell-model to calculate the phonon dispersion curves of 
bulk Sr2Cu02Cl2. In order to ensure dynamical stability 
over the entire BZ, acu-o was increased by %15. 

The shell-model parameters, which include the ionic 
charge Z, the shell charge Y and the intra-ion shell-core 
force constant K, are also listed in Table [I] These pa- 
rameters were adjusted to give a best fit to the four 
measured^i*^ IR-active phonon frequencies at the T- 
point. The agreement between measured and calculated 
frequencies is reasonably good, as can be seen in Table 
ITT1 and the discrepancy between measurement and calcu- 
lation is at most %5. 

The eigenvectors and frequencies obtained from the 
shell- model for the transverse optic (TO) phonon modes 
of Sr2Cu02Cl2 at q = are depicted in Figure EJ The 
four E u modes are the IR-active modes listed in Ta- 
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FIG. 6: The eigenvectors and frequencies of the transverse 
optic (TO) phonon modes of Sr2Cu02Ci2 at q — as obtained 
from the shell- model. All the frequencies are in units of cm -1 . 



blc [n] and are doubly degenerate. The atomic displace- 
ments corresponding to these modes are in agreement 
with mode assignments of Ref. 41. E u (l) is a Cu-0 
stretching mode, E u {2) is a Cu-0 bending mode, E u (3) 
is the translational vibration of a Sr-atom layer against 
the octahedron and E u (4) is the apical bending mode of 
CI against the Cu-0 unil^i^. The non-degenerate A^ u 
modes are also IR-active and correspond to the atomic 
displacements along z-direction. The doubly-degenerate 
Eg and non-degenerate A lg modes are Raman-active and 
£?2m is silent. 

The bulk BZ for body-centered tetragonal structure is 
shown in Figure |7| with the labelling of high-symmetry 
points and directions. Bulk phonon dispersion curves, 
calculated along A, E and TN directions using the shell- 
model of TableHI are shown in Figure[S] The longitudinal 
acoustic (LA) and transverse acoustic (TA) modes along 
all three directions are also labeled. It can be inferred 
form these dispersion curves that the dynamics of the 
Sr2Cu02Cl2 bulk is stable throughout the BZ without 
any imaginary frequency. 



B. Slab Lattice Dynamics 

Our studies of the surface lattice dynamics of 
Sr 2 Cu0 2 Cl2 were also based on the shell- model. In these 
studies we employed slab geometries comprised of SrCl 
- Cu02 - SrCl repeat unit with 21 and 33 layers and 




FIG. 7: The bulk BZ for body-centered tetragonal structure 
with the labelling of symmetry directions. 



with periodic boundary conditions in the plane of the 
slab. The slab was bounded by two SrCl surface layers. 
The total number of degrees of freedom in the 33-layer 
slab is 231. The projection of the bulk BZ on the sur- 
face is illustrated in Figure[!|| Projections of all the bulk 
phonon modes that lie inside the SBZ contribute to the 
bulk bands in slab dispersion curves. 

Figure shows the results of the slab calculations 
along both (100) and (110) directions of the SBZ where 
unmodified bulk shell-model parameters (Table [I| were 
used throughout the slab. The grey areas and curves in 
Figures 110b and 110b are projections of the bulk modes 
on the SBZ (bulk bands). In Figure ITUk solid black lines 
identify dispersion curves for shear horizontal (SH) sur- 
face phonon polarizations. These modes are polarized 
perpendicular to the sagittal plane (defined by the sur- 
face normal and the phonon wave vector) and parallel 
to the surface; they are odd with respect to reflection in 
the sagittal plane. Because in our scattering geometry 
the sagittal and scattering planes coincide and contain a 
high symmetry surface direction, SH modes would not be 
detected. In Figure HUb. solid black lines correspond to 
surface modes which are polarized in the sagittal plane 
(SP modes) and the experimental points are shown by 
filled squares. 

From the surface dispersion curves of Figure 1101 it is 
clear that all the SrCl surface modes are below 200 cm -1 . 
The high energy bulk bands, labeled as longitudinal optic 
(LO) and transverse optic (TO) bands in FigurelTUk. cor- 
respond to the vibrations of Cu and O atoms in the bulk 
layers, consistent with the results from bulk dispersion 
curves. 

Below 200 cm -1 , the agreement between the experi- 
mental points and the surface modes is reasonably good. 
In both directions there are experimental points which 
clearly fall along the surface SP modes. One example 
is the SP mode around 115 cm -1 along (100) (labeled 
by (1) in Figure HUb ). which corresponds to vibration of 
surface Sr atom along z direction and CI atom along x 
direction. Another example is the mode labeled by (2) 
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FIG. 8: Phonon dispersion curves along three symmetry directions in the bulk BZ calculated by the shell-model of Table [I] 
The four measured—*— IR-active E u frequencies at V point are labeled by asterisks. 







r z 


' A 




<100> 



(a) (b) 

FIG. 9: (a) SBZ and (b) projection of the bulk BZ on the 
surface (dashed-lines) . 



along (100) direction which corresponds to the displace- 
ment of the surface CI atom along z. Another area where 
there is agreement between theory and experiment is the 
region with a high density of bulk bands from ~ 30 to 
70 cm -1 , close to the T point along (100). Because of 
the high density of bulk states in this region, a large 
number of experimental points have also been detected. 
Along (110) direction, mode (3) is associated with the 
vibrations of the surface Sr atom along z and mode (4) 
corresponds to displacements of both Sr and CI surface 
atoms in the xy plane. 

Along both directions, the acoustic Raylcigh modes are 
labeled by R. The scarcity of experimental points along 
the Rayleigh modes is attributed to the limitation of the 
energy resolution that has prevented modes with energies 
lower than ~ 25 cm -1 to be clearly identified. 



Because of the systematic agreement between the SrCl 
surface modes and experimental dispersion points, no fur- 
ther adjustments of model parameters were necessary. 

In order to confirm that the cleaved (001) surface of 
Sr2Cu02Ci2 consists of only SrCl layers and not Cu02, 
we implemented the slab calculations for a 21-layer slab 
with a CuC>2 surface termination. In order to maintain 
the charge neutrality, the slab had to be made asym- 
metric, i.e. the other surface of the slab was a SrCl 
layer. Employing bulk shell-model parameters resulted 
in instabilities of several modes away from T. Efforts to 
completely suppress these instabilities were not success- 
ful. However, we found that setting &cu-o m the surface 
layer to 5.2 A -1 resulted in confining the instabilities to 
the lowest lying modes. In these calculations the optic 
modes of the Cu02 surface were stable as shown in Fig- 
ure ^2 for SP surface modes along (100) direction. The 
instability of the lowest lying acoustic modes away from 
the f point is clear in Figure ITTb. 

By contrast to the results of slab calculation for 
SrCl surface, we did not find systematic agreement be- 
tween the experimental data and calculated Cu02 surface 
modes. For example no experimental points exist close to 
the surface modes labeled (1), (2) and (3) in Figure ITTI 
Mode (1) corresponds to the vibration of 0(2) surface 
atoms in the x direction, mode (2) corresponds to the 
displacements of the surface 0(1) and 0(2) atoms along 
z and mode (3) is associated with vibrations of surface 
Cu and 0(2) atoms along z and 0(1) atoms in the xz 
plane. The coincidental agreement of a few points with 
some Cu02 surface modes (e.g. mode (4) around 125 
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FIG. 10: Results of the slab calculation for a 33-layer slab of Sr2Cu02Cl2 with SrCl surface termination along (100) and (110) 
directions. The grey curves show the bulk bands. In (a) black solid lines correspond to surface modes with SH polarization, 
while in (b) the black lines show SP surface modes. The filled squares are the experimental data obtained by inelastic He 
scattering. Modes labeled by (1), (2), (3) and (4) are discussed in the text. The acoustic Rayleigh modes are labeled by R. 



cm 1 ), is more accidental than systematic. 

V. CONCLUSION 

Cleaved (001) surfaces of the layered perovskite 
S^CuO^Ck" were prepared in UHV, and the structure 
and phonon dispersions of the surface were studied ex- 
perimentally using elastic and inelastic HAS along high- 
symmetry directions of the SBZ. The position of the 
Bragg peaks in the diffraction patterns, obtained for dif- 
ferent incident angles, revealed a surface periodicity con- 
sistent with bulk termination, without any evidence of 
surface reconstruction. This fact, along with previous 
LEED and x-ray photoemission results2^Si2& confirms 
that the resulting surface is non-polar and stable which 
favors a SrCl surface termination. A corrugation func- 
tion for the SrCl surface was constructed with the aid of 
eikonal approximation using measured diffraction peak 
intensities and was in agreement with a corrugation func- 
tion obtained from surface charge density calculations. 

Bulk lattice dynamical calculations based on the shell- 
model were carried out by transferring some of the poten- 
tial parameters from existing models for other cuprates 
like La2Cu04 and obtaining the rest by solving static 
equilibrium equations. The shell-model parameters were 



adjusted to fit the four experimentally measured IR- 
active phonon frequencies at q = 0. No instability were 
observed over the entire BZ. 

Surface dynamics was investigated using slab geometry 
and the same shell-model. The calculations for a 33-layer 
slab with SrCl surfaces and a 21-layer slab with a CuC>2 
surface revealed two important features. First that the 
experimental surface phonon dispersions obtained by in- 
elastic HAS agreed quite reasonably and systematically 
with a SrCl terminated slab model. Second, the shell- 
model for a Cu02 terminated slab could not be stabilized, 
despite all efforts. This and the fact that no experimen- 
tal dispersion points lie systematically along the Cu02 
surface modes is another sign that the surface layer is 
SrCl and Cu02 sheets do not exist on the surface. 
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FIG. 11: Results of the slab calculation for a 21-layer slab of Sr2Cu02Cl2 with CUO2 surface termination along (100) direction. 
The grey curves are bulk bands, the black solid lines are the surface modes with the SP symmetry and the filled squares are 
the experimental data. Modes labeled by (1), (2), (3) and (4) are discussed in the text, (b) is the expansion of (a) from 75 to 
225 cm -1 . 
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